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We present an 8.4-GHZ radio imaging study of an optically selected sample of early type Seyfert 
galaxies; comparison of images taken at two epochs reveals possible variation in the nuclear radio 
flux density in five of them over a seven year period. It is shown that there is a possible correlation 
between the presence of nuclear radio variability and the absence of hundred parsec-scale radio 
emission, analagous with radio-loud AGN. Our results suggest that all Seyferts may exhibit vari- 
ation in their nuclear radio flux density at 8.4 GHz, but that variability is more easily recognized 
in compact sources in which emission from the variable nucleus is not diluted by unresolved, 
constant flux density radio jet emission within the central 50 pc. Taken in combination with 
other Seyfert properties, these results suggest a paradigm of intermittent periods of quiescence 
and nuclear outburst across the Seyfert population and demonstrate the importance of investigat- 
ing the temporal domain at radio wavelengths, which remains completely unexplored for faint 
radio-quiet AGN. Discovery of intermittent activity and radio flares has important implications 
for the AGN duty cycles and the origin of Ultra-High Energy Cosmic Rays. New radio facilities, 
such as the EVLA/VLBA, eMERLIN, LOFAR and eVLBI on the EVN, will revolutionise the 
study of radio-quiet AGN; in particular, the combination of increased sensitivity and sampling 
rate with high-angular resolution and automatic data reduction will open up the transient radio 
sky and bring major future breakthroughs. 
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1. Introduction 

Variability of nuclear flux density at all wavelengths across the electromagnetic spectrum has 
long been recognised as a defining characteristic of Active Galactic Nuclei (AGN) |jljj. The short 
variability timescales measured at X-ray, optical and UV wavelengths, ranging from days to years 
depending on observing wavelength and intrinsic AGN luminosity, support the standard AGN 
paradigm of nuclear emission originating from compact regions powered by accretion of matter 
by a central supermassive black hole, rather than standard stellar processes BJ|. 

Variability at radio wavelengths is most marked in powerful radio-loud AGN. The largest, 
most rapid variations in nuclear non-thermal continuum emission are seen in BL Lacs and core- 
dominated quasars These variations are explained primarily by Doppler boosting of the nuclear 
emission by highly-relativistic jets viewed at small angles to the line of sight and are consistent with 
radio-loud Unification Schemes [Q]. Flares in the radio light curves might correspond to ejection 
of new relativistic components or emergence of shocks in the underlying flow 

Seyfert nuclei, classified as radio quiet though not radio silent, have been known for some 
time to show variability in their continuum and line emission at X-ray, optical and UV wavelengths 
[0,|8|], providing constraints on nuclear absorbing column variations, photon reprocessing, structure 
and dynamics of the Broad Line Region and mass of the central object. However, although ten 
times more common than their radio-loud counterparts, Seyferts are 100 to 1000 times weaker at 
radio wavelengths [^] ; therefore, little is known about the radio variability of their nuclei as only 
a small number of 'interesting' Seyferts have been systematically monitored at radio frequencies, 



with some studies motivated by serendipitous discovery of radio flares [|10|]. 

Some Seyferts have small radio jets, < 1 kpc in extent [O, [l2|] ; with a constant flux density and 
steep spectral index, this extended radio emission may obscure fainter flat-spectrum emission asso- 



ciated with the AGN itself unless observed at sub-arcsec angular resolutions [ 13, 14]. Identification 
of compact flat-spectrum radio nuclei with high brightness temperatures using the VLBA, showed 
unambiguously that Seyfert nuclei are black-hole driven central engines - not compact nuclear star- 



bursts Q15Q - and, coupled with the identification of a highly collimated radio jet in NGC 4151 (Fig. 



HI) [O, 16], suggests that radio variability in Seyferts may be common. 



2. A VLA Search for Nuclear Radio Variability In Seyferts 

We conducted a two-epoch study of 8.4-GHz nuclear radio emission in a small, but complete 
sample of 1 1 nearby optically-selected Seyferts in which the nuclear radio emission is imaged 
at angular resolutions of ~0.2" (<50 pc); the original goal of the study was to obtain high quality 
radio images to investigate the impact of Seyfert radio jets on the interstellar medium by comparing 
extended radio structures with extended ionised gas distributions and excitations in the Narrow Line 
Region as inferred from images with the Hubble Space Telescope (HST) [|l7|]. This comparison is 



presented in Mundell et al. (2009) [18], including new radio images of MCG -5-23-16. However, 
by more reliably separating nuclear and extended radio components and comparing with previous 
nuclear flux density measurements we serendipitously discovered a fraction of Seyfert nuclei with 



variable nuclear radio flux densities [18] 
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3. Radio Variability, Seyfert Type and Radio Morphology 

We selected all 12 Seyfert galaxies in early type hosts with my < 14.5, cz < 3000 km s _1 and 
declination suitable for VLA observations from Mulchaey's sample of early-type Seyfert galaxies 



[19]. The optical continuum and ionized gas properties were studied with the HST [1170; all ob 



jects, except MCG -5-23-16, were observed with the VLA in 1993, at 1.4 and 8.4 GHz, as part 



of the larger VLA snapshot study of [ |20| | (hereafter N99) providing a comparison epoch for our 
observations, which were conducted using the VLA in A configuration at 8.4 GHz on 1999 August 
20/21. Full details of data reduction and analysis - including careful and important corrections for 
tropospheric phase decorrelation - are given in Mundell et al. 2009 [O]. Data presented in N99 



were reprocessed using the method presented in [18] to allow rigorous comparison of flux density 



measurements of the weakest sources. The results are as follows: 



• All twelve Seyferts are detected with the VLA; five nuclei show a variation in their nuclear 
flux density between 1992 and 1999 ( NGC 2110, 3081, 5273 and MCG -6-30-15 declined 
and NGC 4117 increased). The remaining six - Mrk 607, 620, NGC 1386, 3516, 4968 and 
7465, remained constant (one was not observed in 1992). 

• No correlation is found between the detection of variability and core luminosity, which lies 
in the range 8.2 x 10 18 < L(core) 84 ghz < 5.9 x 10 21 W Hz -1 for this sample. 

• No correlation is found between Seyfert type and the detection of radio variability. Instead, 
the compactness of the radio emission seems to be a good indicator of variability. All variable 
nuclei, apart from NGC 21 10, are compact and at best marginally resolved, while nuclei with 
associated extended radio emission have remained constant. 

• NGC 2110 (see Fig. is the only Seyfert in this sample with significant extended radio 
emission from well-defined radio jets and a variable nucleus, suggesting that the intrinsic 
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Figure 2: VLA A-configuration 8.4-GHz radio image of NGC 2110 taken from Mundell et al. (2009). 

nuclear variability is higher than the observed ~38%. The observed nuclear flux variability 
indicates significant nuclear structural changes on scales smaller than the VLA beamsize (i.e. 
within the central 15 pc), between the two epochs, possibly due to the appearance and fading 
of new components or shocks in the jet. This is consistent with sub-parsec scale structures 
identified by Jl5| ] that consist of a high brightness temperature (> 6.0 x 10 K), synchrotron 
self-absorbed point source, extended jet emission 0. 1 pc north and south of the nucleus, and 
a discrete component ~0.3 pc north of the nucleus, possibly the first knot in the northern jet. 

• The serendipitous discovery of variability in ~50% of our sample, which was not selected 
on any given radio property or known variability, suggests that radio variability in Seyfert 
nuclei is a common phenomenon and the measured amplitudes of variation are likely to be 
lower limits to peak-to-peak variations within this period. This is consistent with studies 
of III Zw 2, which flares approximately every five years [pT|] - a timescale consistent with 
the variability detected for our sample, given that two-epoch observations are more likely to 
detect fading sources if the decay timescale is longer than the brightening timescale. 



All Seyferts might exhibit variations in their nuclear radio flux density at 8.4 GHz, with it 
being more easily recognised in compact sources in which emission from the variable nucleus 
is not mixed with unresolved, constant flux density radio-jet emission within the ~0.2" VLA 
beam ( central < 50 pc for V. svs <4000 km s _1 ). This is consistent with an increased detection 
rate of flat spectrum nuclear radio emission in Seyferts imaged at VLBI resolution [15, E3f|, 



further emphasising the importance of isolating nuclear emission from jet emission. Linear 
resolution limitations for studies of more distant AGN such as radio quiet quasars (RQQ) 
are more accute. Barvainis et al. p2[ ] find few differences between RLQ and RQQ other 
than an unexpectedly weak dependence of radio variability on nuclear spectral index, with 
RQQs showing steeper spectral indices than RLQ. If radio quiet quasars have extended, but 
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unresolved steep spectrum emission similar to Seyferts, their nuclear variability and spectral 
index might have been underestimated, thus increasing their similarity to RL AGN. 

• If flares in radio light curves correspond to ejection of new relativistic components or emer- 
gence of shocks in the underlying flow, sensitive systematic monitoring of a larger sample, 
combined with subarcsecond resolution imaging, may confirm that Seyferts - as black-hole 
driven AGN - have the capacity to accelerate relativistic jets during radio flares despite hav- 
ing radio jets that are intrinsically non-relativistic during quiescence. 



4. Radio Flares and Relativistic Jets 



Flux variability and apparent superluminal motion are defining characteristics of bulk rela- 



tivistic motion in extragalactic jets. Long-term monitoring of radio-loud AGN [ |24| ] demonstrated 



the existence of stochastic, non-random flare events at cm/mm wavelengths, with mean timescales 



2 years for BL Lacs and QSOs [24]. In standard jet-shock models, increases in total flux density 



are attributed to the passage of shocks on pc scales in collimated, relativistic flows and, in some 



cases, radio flaring corresponds to propagation of individual components in VLBI maps [25]. 

In contrast, debate continues over the nature of ejected radio plasma in Seyferts; the presence 
of bulk relativistic motion has not yet been demonstrated in Seyfert jets, with existing proper motion 
studies measuring apparent radio component speeds V app < 0.25 c [[!(]]. The primary limiting factor 
in establishing the radio variability characteristics of Seyferts and their radio jets has been the lack 
of high resolution, long-term monitoring surveys of large samples of Seyferts, whose radio nuclei 



are faint in quiescence so require careful, phase-referenced observations (e.g. [|12J). 

Despite the small number of Seyferts studied in detail so far, there is growing evidence that 
their nuclear activity is characterised by intermittent quiescent phases punctuated by outbursts [^6|]. 
Serendipitous discovery of radio flares in Mrk 348 and III Zw 2 - suggest that Seyferts can accel- 



erate their jets to large speeds during an outburst [27 ]. Most notable was a flare in III Zw 2 in 1998 



that showed a ~30-fold flux increase within 2 years and a highly-inverted radio spectrum peaking 
at millimeter wavelengths; similar flares occurred at earlier epochs over the last ~25 years at lower 
frequencies [p7|]. A sudden change in spectral shape in 1999 indicated a structural change in the 
nucleus - confirmed by VLBI imaging of the nucleus that revealed a newly ejected radio component 
expanding at an apparent superluminal speed V app > 1.25c [Q], in stark contrast to the upper limit 
of 0.04c measured during quiescence. Discovery of a radio flare in NGC 2110 and the variability 
of compact Seyfert nuclei in our sample suggests this may be a common phenmenon. 

4.1 Radio Flares in Seyferts as the Origin of Ultra High-Energy Cosmic Rays 

The Auger Collaboration recently reported a significant anisotropy in the arrival directions of 
Ultra-High-Energy Cosmic Rays (UHECRs) and a possible correlation of events with catalogued 
local AGN; this is consistent with the theoretical prediction that energy losses due to photopion 
production on the cosmic microwave background limit the sources of the highest energy CRs to 
D<100 Mpc. Although particle acceleration in jets and lobes of local radio galaxies offer promis- 
ing candidates for some of the CR events, many remain controversial and debate centres around 
whether the coincidence between the relatively large number of local AGN compared with the small 
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number of CR events is mere chance alignment [28]. Theoretically, very modest, but intermittent 



jet activity in low luminosity AGN, such as Seyferts, could provide the acceleration mechanism 



for UHECRs [ |29| , |30p. In particular, acceleration of UHECRs and the survival of energetic heavy 
nuclei are possible in the parsec-scale weak jets that are typically observed in radio-quiet AGN. 

Routine, high quality, high angular resolution imaging of a well-selected sample of Seyferts 
with current and future radio facilities, such as EVLA, eMERLIN, VLBA and LOFAR coupled with 
multi-messenger follow-ups, would provide valuable input to the debate on the temporal properties 
of Seyferts, as well as investigating the wider possibility that duty cycles and jet production in local 
AGN may be associated with UHECRs. 
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